Setup for the Josephson effect 



In this section we describe in greater detail the setup we considered in the linear response calculation, relaxing 
some of the conditions assumed in the article. Altough we use the notation and terminology of Figure la of the main 
article the description is equivalent for example with the system of Figure lb of the main article. 

In order to accommodate a supcrfluid with balanced pairing, we require a balanced number of particles within 
the supcrfluid states i.e. Ni = N2 and N3 = N4. We further assume that for any tunneling in the system the 
relation SNi/Ni <C 1 holds, so that we do not have to take into account the internal dynamics of the supcrfluids. 
In particular, this assumption implies that the supcrfluid states are not destroyed because of the transitions. The 
tunneling link between the hyperfme states is provided by two radio frequency (RF) fields Ei3(r, t) and E24(r,i). 
The field Ey(r, t) drives transitions between the states \i) and \j). For simplicity we take these fields as a single 
mode with frequency Vij and wave vector ky . The strength of the coupling between the electric field Ejj (r, t) and 
the hyperfme transition between states \i) and \j) is given by the Rabi frequency 



Ey(r), 



CD 



where djj is the electric dipole moment of the transition. Throughout the treatment we assume that the RF fields 
are near resonance and that the intensity of the radiation is low. With these assumptions we may use the rotating 
wave approximation [T]. As a result, it is possible to transform the time dependence of the field in the rotating frame 
into a shift in the transition frequency Sij, called detuning, defined as 

Sij = vij — ujij, (2) 

where is the frequency of the transition. 

Within the same supcrfluid state the hyperfme states will naturally have the same chemical potential, therefore 
A'i = P"2 and /13 = (14,. In the article we further assumed fii = /i3, but here we remove this constraint. Including this 
difference in the chemical potential between the two superfluids into the RF detuning, we get the modified detuning 
of the hyperfme transition 

Sij = + fj,j - Hi. (3) 

This quantity describes the energy gain or loss associated with a particular hyperfme transition without many-body 
effects. 
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Figure 1: The energy levels of the system. In the unperturbed system the chemical potentials are equal within the 
BCS state: /ii = ^2 and ^3 = (14. The detunings ^13 and S24 of the RF fields can create an effective asymmetry in 
the energies. 



The energy levels of the system in question are shown in Figure [T] Here we make use of the assumption that the 
applied RF field does not alter the internal dynamics of the superfluids so that we may split the detuning Sij between 
the states \i) and \j) freely. We choose a symmetric splitting, although we will eventually see that this choice does 
not affect the results. 



Linear response description 

The system is described by the Hamiltonian 

H = Hq + Hrf . (4) 



Here we have defined Hq as 
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Ho=j dr£>!(r) (-1 

+ \f rfr^^i^(r)V](r)^(r)Vi(r). (5) 



The external perturbation is given by 
where we have 



H RF = H C + H d , (6) 



B - = / dr E fi ^J( r )^(r) (7) 
ff d = [ dT^WrtHWii*)* ( g ) 



with the definitions Wi = fii + 5is/2, W3 = ^3 — J13/2 and similarly for W2 and W4. The inclusion of the chemical 
potential in these expressions is just a mathematical trick to aid the forthcoming calculation. For RF fields it is a 
very good approximation to take the Rabi frequency Q as independent of position. 

The operator for current density to one hyperfine state is given by the time derivative of the particle density 
operator, which is obtained using the Heisenberg equation of motion 

j t N i (v 1 )=i[H,N i (v 1 )]. (9) 

Only H c fails to commute with the number operator. Hence the Heisenberg equation of motion for the number 
operator becomes 

±N i (r 1 )=i[H c (t),N i (r 1 )) 

= i E/ dv^^^^Aivi)-^^!)^))- (10) 
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The physical current density h(yi) is the thermodynamic average of the current density operator defined as Ij(rii) = 
/TVj(ri)^, which we then need to calculate. However, the exact calculation of this quantity is not feasible, since time 

evolution of the operators in the Heisenberg picture depends on the coupling H c . Hence, we need to resort to 
an approximate solution with respect to the coupling strength. At this point it is most convenient to make a 
transformation to the interaction picture with respect to H c . For a general operator A(t) this transformation is given 
by 

A(t)=S(t,t )- 1 A(t )S(t,t ), (11) 

in which to gives the initial condition and 

S(t,t Q ) = Tcxp ^-i J dt' H c {t')j . (12) 

Assuming that H c is a weak perturbation, it is a good approximation to linearise this transformation with respect 
to H c , in which case we have 

t 

S(t,t ) = 1-i I dt'H c (t'), (13) 



and consequently A(t) is linearised to 



A(t) = A(t )-i J dt' [A(t),H c (t')]. 

to 

The thermodynamic expectation value of a generic operator A(t) then follows the Kubo formula 



(14) 



(A(t)) = (A(t ))-i J dt' ([A(t),H c (l/)}), 

in 



(15) 



Here to is the reference time for the initial condition. Since we assume a system initially in equilibrium, there is no 
current at the time to. Therefore, upon replacing A with the current operator N, the first term in the expression 
(|15l) vanishes. In order to obtain the Kubo formula in the frequency domain, it is a common strategy to take the 
limit to ~oo for the initial condition. 

Inserting equation (fTCTj) into equation (TT5)) then yields the current density 

oo 

/i(rit) = i f dt'J2 [drzdrsdrJeMiW-WJt + i (Wj - W k ) t'} 
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where the retarded linear response function Lij k i(Y\t, r 2 i, r 3 i', r 4 i') is defined as 

L ljk i(v 1 t,v 2 t,v 3 t',r 4 t') = -^(t-t / )([V'feM)V ) i(rit),V'i t (r4t')V' J '(r 2 *')]) 



(16) 
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In equation (|16|) we have used standard algebraic manipulations to separate out the effect of the quantities Wi from 
the remaining time evolution of the field operators. For the hyperfine species |1), equation (|16|) can be rewritten as 



J 13 (r 1 t) = /f 3 (r 1 t)+/ 1 J 3 (r 1 t), 
where the standard single particle (quasiparticle) current Jf 3 is defined as 
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and the Josephson current J/ 3 as 



I( z (v x t) = -2 / rfr 2 Im ^e-^+^n^il^ 

x Li234(ri - r 2 , -5 24 + ii^ 
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In Eqs. (fl9l [20]) we have identified the Fourier transformation in time for L. This was the reason for separating the 
quantities Wi from the time-evolution. In order to quarantee the convergence of the transformation, an infinitesimal 
(positive) convergence factor r) + has been introduced. Thus, we see that the modified detunings defined already in 
equation (j3J) now appear in the expression for the current. Here we may identify the spatial integration as a Fourier 
transformation with respect to zero momentum. Furthermore, in a homogeneous system we may drop the variable 
i"i since the current is not a function of position. Therefore we obtain 



7f 3 (t) = -21m [^ 3 f> 13 Li33i(p = 0A 3 +iv + )\ , 
I( 3 {t) = -21m le-^+^ni^L^ip = 0, -<S 24 + 



(21) 



We note here that the definitions of single-particle and Josephson currents agree with their definitions in the con- 
ventional case of two superconductors when J13 = 5 2 4- 



At this point in our calculation we have isolated the effect of the perturbation on the initial system in the linear 
response approximation. Hereafter we only need to work with linear response functions which depend only on the 
properties of the unperturbed system. 

Kadanoff-Baym formalism 

Our goal here is the calculation of the retarded linear response function L. To this aim we will resort to the 
Kadanoff-Baym formalism [2} 13] to which we refer also as the self-consistent method. The interest of this method for 
our purposes lies in the fact that it allows the inclusion of the cross-interaction effects in the transport properties, 
and that, on more general grounds, the linear response function L will obey the same conservation laws that are 
obeyed by the single-particle Green's function Q. 

In the Kadanoff-Baym formalism, when the general case of a time and position dependent external perturbation 
is considered, L is obtained from Q by functional differentiation with respect to the external field 
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Here we are working in complex time Matsubara formalism, which is a standard technique for dealing with finite 
temperature Green's functions. The real-time retarded linear response function can be obtained from L in the 
Matsubara formalism by analytical continuation. 

For the sake of brevity, we truncate the variables such as r,; and the spin-index in the following two equations 
to just the index i. Furthermore, we denote the variables of integration with a bar over the variable. The essential 
starting point is the equation of motion for the single particle Green's function 



^(l, 1)0(1,1') = 5(1-1') + J fi(l, 1)0(1,1') + I S(l, 1)0(1,1'), 



(23) 



Here Qq 1 is the inverse free propagator, which accounts for the kinetic part of the equation of motion, and the self- 
energy £ contains all the interaction effects in the system. On the level of principle one should solve this equation 
for an arbitrary perturbation 51 to obtain Q and then further calculate L. This is in practice unfeasible, and in the 
Kadanoff-Baym method one circumvents this problem by taking the variational derivative of the equation of motion 
at O = 0. This leads to an implicit equation for L 



L(12,l'2')= J 0(l,3) n =o0(4,l')n=o 
+ / 0(l,3)o = o0(4,l')o= o 
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which is easier to solve. In fact, in a homogeneous system we may utilize Fourier transformations leading to a matrix 
equation for L in momentum space. 

In the following calculation we consider the case where U13 ^ 0; there is no conceptual difference when also other 
cross interactions are considered different from zero. If the BCS approximation is considered for the unperturbed 
system, the explicit form of the linear response function becomes 
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for the current between states |2) and |4). 



The terms II in the previous equations are 

n<?(p,w) = -^^<m(q,x)033(q- p^x - w), 
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Here J- is the so called anomalous Green's function, which expresses the pair correlation of the supcrfluid. In the 
main article we study a case where the two supcrfluids are identical, in which case we have Tig = Tig. 

The expression for the linear response function leads to the following results for the single-particle and Joscphson 
currents 
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The phase factors ^(^13) and tp(S 24 ) which appear in equations (|32[) and (|33|) are the initial complex phase of the 
retarded linear response function L. In the article we have left out the minus sign from the detunings in the terms II 
since in the case of identical supcrfluids, which we consider in the article, the critical current becomes a symmetric 
function of the detuning. 



Identical superfluids 

Here we discuss the results obtained in the case of identical supcrfluids. In this case the notation simplifies 
as £13 = <5i3 + /Lt 3 — fix = <5i3 and 5 2 ± = 5 24 + M4 — ^2 = ^24 since the chemical potentials are equal. From the 
Kadanoff-Baym formalism, in the case where cross-interactions are zero, we have 

/faM) = -2|^i 3 | 2 Im [n e (p = 0, <5i 3 + irj + )] , 
iL(rit) = -2\n 2i \ 2 Im [n g (p = 0, S 2i + , 

and 

I&(r x t) - 2 |Sl 13 f7 24 n^(p = 0, -S 2i + ir) + )\ sin ((<S 13 + 5 2i )t + tp^fa)) , 
Ii(r x t) = 2 |f> 13 f> 24 IV(p = 0, -* 13 + sin((<5 13 + S 24 )t + <p2(5 13 )) . (35) 

The single-particle current vanishes when \Sij\ < 2A 12 (Fig. [5]). In this range only the Josephson current exists, 
as it is intuitively clear, since the single-particle current requires the breaking of Cooper pairs. At |<5y| = 2A 12 the 

detuning equals the minimum of the transition energy Ei 2 (q) + E^q) , where E a p(q) = J (\q\ 2 /2m — fi a ) 2 + \A a p\ 2 . 
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This occurs at the point -2— = /i i.e. at the Fermi level, where the BCS density of states is divergent, leading to 
a divergence in the dependence of the currents on the detuning. In the case of the critical Josephson current, this 
divergence is commonly known as the Riedel peak. 

To illustrate the spin-asymmetry further, consider the results of Fig. [2] Notice that we plot only Jf^ since I 24 
is the same function, the difference being only that ![% depends on 5 24 whereas I 24 depends on <5i3 . For instance 
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Figure 2: (left to right). The single-particle current /-f 3 and the critical Josephson current /f^ vs. S13 and S24 for 
identical BCS states. U12 = U34 = — 2.0Ep/n. T is in the units of Fermi temperature Tp and in the temperature 
range of the figure, the gap is A ~ OAEp. 
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Figure 3: (left to right). The single-particle current 2f 3 and the critical Josephson current vs. S13 and S24 for 
unidentical BCS states. U12 = — l.OEp/n, U34 = —S.OEp/n, T is in the units of Fermi temperature Tp. 



taking S 13 /E F = 0.4 and S 24 /E F = 0.5 at T/T F = 0.08 then leads to lf 3 {524) / 1? 4 (5 13 ) = 1.14, in other words, a 14 
% difference in the amplitude of the Josephson current for each spin component. 

Finite temperature and cross interactions 

It is instructive and of certain interest, from the experimental point of view, to discuss how the results presented 
in the article arc affected by temperature and different interaction strengths. In order to properly observe finite 
temperature effects, we consider the slightly more complicated case of two different BCS states. Experimentally, the 
condensates 1-2 and 3-4 may have different interaction strenghts. We choose U12 — —l.OEp/n and U34 = —h.QEp/n. 
This leads to gaps A12 ~ 0.2Ep and A34 ~ 0.5Ep, which corresponds roughly to systems that have been realised 
experimentally. Also the choice of temperature we use, T = 0.06 . . . O.lTp, is motivated by the lowest temperatures 
reached in the experiments. The results are plotted in Figs. [3] and 2] First of all we notice that all the features 
previously found still dominate the result, even though there are some clear changes. 

We see that for finite T and U12 ^ U 3 4 there can be a finite single particle current even at small detunings. This 
current arises due to thermal excitations in the BCS states. Intuitively, this contribution is generated by Cooper 
pairs that have been already broken because of thermal energy. To be more precise, it is the fact that one of the BCS 
states has more thermal excitations for a given energy than the other that leads to the emergence of the thermal 
single particle current. We did indeed have thermal excitations in the results of the previous section as well, but 
their net contribution summed up to zero. The new peaks in the results (see Figures [3] and [4]) are associated with 
the the contribution of thermal excitations diverging at a particular detuning. The explanation for these divergences 
is similar as for the divergences encountered in the previous section. Here the energy of the tunneling process is 

E\2 (q) — E 3 4(q) (again with E a p(q) = J (^rl^l 2 — M«) 2 + |A Q/ 3| 2 ) and the resonance occurs when the detuning 
coincides with the extremum of this energy. The density of states argument of the previous section holds then here 
as well. Moreover, the thermal single particle current vanishes above this maximum because at larger detunings 
there are no thermal transitions that would conserve both energy and momentum. 
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Figure 4: Single-particle current /-f 3 , and critical Josephson current 1^ vs. Si 3 and S24 for unidcntical BCS states. 
U12 = — l.OEp/n, C/34 = —b.OEp/n, T is in the units of Fermi temperature Tp. 
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Figure 5: Energy levels for the four-state model. 

The perturbative calculation for the four-state system 

The Hamiltonian 

For the perturbative calculation of the four-state system of Figure El we consider the following Hamiltonian 

H = H + H RF , (36) 
where Hq contains the energies of the states and Hrf is the RF coupling. They have the form 
H a = -U {nm 2 + n 3 ni) , 
H RF = Qi 3 c\c 3 + h.c. + Q, 2i c\c A + h.c. + ^ (m - n 3 ) + (n 2 - "-4) , (37) 

where U is the energy difference between the paired and unpaired states. We have kept the notation here as similar 
as possible to the initial description in the article. Obviously one could reformulate this four-state system in a more 
compact form. 

Perturbative calculation 

Now let us study our model system analytically using perturbation theory. The dynamics of the system can 
be described by four states: the paired states \<f>\) = |12), \4>\\) = 1 34) - corresponding to the left and right well 
respectively in the double well description - and the states |14), |23) (see Figured]). 



In order to find the lowest order contribution of the RF-couplings we have to employ the time-dependent pertur- 
bation theory to the second order in the RF-couplings. As an initial state we consider an arbitrary superposition of 
the paired states 

\<h) =ao\<h) +A)|<M- (38) 

We point out already that in the ao = or (3q = case no Josephson current will appear. Calculating the occupation 
number Ni of the hyperfine state in the state \<j>(t)) = exp [—iHt] |</>q) gives 

(0(t)|JVi|0(t)> = \a \ 2 (MVWiMt)) + \Po\ 2 (Mt)Wi\Mt)) 

+ao(3* (Mt)\Ni\Mt)) + a*oPa (<M*)l^i|M*)) ■ ( 39 ) 

The perturbative expansion for |</>i(f)) and \4>u(t)) reads 

\Mt)) = |12> + j$(t) |14) + 8$(t) |23) + [a$ 12 (t) + ag\ 2 (t)} 1 12) + [0&> 12 (f) + (3g\ 2 (t)] |34) 
\<ki(t)) = 134) + yg\t) |14) + 8<£(t) |23) + [a^ M (t) + a% M (t)] |12) + [/3$ 34 (i) + pg{ 4 (t)} |34) (40) 

in equation (|40[) the superscript represents the perturbative order, the subscript the path corresponding to the process 
considered and the Greek letter its starting state (with the convention a <H> |12), (3 <H> |34),7 <-> |14),<5 o |23)). For 

(2) 

instance (3\ 412 corresponds to the second-order process 1 34) — > 1 14) — >■ |12). The coefficients that we need for 
calculating (JVi 
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are obtained within time-dependent perturbation theory, considering Hq as the initial state and switching on Hrf 
at t = 0. In the equation above we have defined 
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Since we are interested in the time evolution of the expectation value of the particle number in state |1) up to 

(2) (2) 

second order in fi, the only relevant second-order coefficients are 0:^4 34 (f) and 34(f)- This is due to the fact 
that the states corresponding to other 2 nd order coefficients do not contribute to (Ni) through the scalar product 
with 1 12) (0 th order term). The calculation of the terms (</>i(f)|ATi|0i(f)) and (</>n(f)|iVi|</>ii(f)) give the quasiparticlc 
current (Ifj in equation (1) of the main article) when the proper normalization is taken into account. 

In the term ao/?o (0n(f)l-^i|<Mf)) + a oA) (<fri(t)\Ni\<fru(t)) we have contributions both from the scalar product 
between 1 st order terms ((14 |. . . 14)) 
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and from the scalar product between the Q th -2 nd order terms ((12 |. . . 12)) 
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Here ip stands for the initial phase difference between ao and (3o i.e. ccqPo — aoA) cxp(iip). In Equations (|4T>)) 
and (|47[) . the components oscillating with frequency (W1434 — uj^ i 2 ) correspond to the contributions giving rise to 
the Josephson current, while components oscillating at W1434, W2312 represent a contribution to the quasi-particle 
current, and the terms in W1412 cancel out. 

If we now focus on the component of (Ni) oscillating at the Josephson frequency and replace y with their 
values in terms of interaction energies and detunings as given in equations (|42p and (f4"3")l. we obtain 



(Nl)j = 2^13^24 
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cos[(5 13 + S 2 i)t + cp}. (48) 



_(U + S 13 )(S 13 + S 24 ) (U + S 24 )(S 13 + 5 2A ) (U - S 24 )(U + S 13 ) 
Note that we would also have a contribution to the Josephson current coming from the scalar product between the 

C2) (21 (2) (2) 

terms containing a\^ 12 (t) + cx 2 3 12(f) an d Q: i4 34 W + a 2s 34(^)1 associated with two different paths, via two different 
internal states, for the tunneling of Cooper pairs, but this is only a fourth order process. Such processes are present 
even if one of the paired states would initially be empty 

Differentiating Equation (|4"8"|) with respect to time, we have 



//g = 20 13 n 24 |ao/3o| [M pair + M sing i e ] sin[(<5i 3 + S 24 )t + <p], (49) 

where M pa j r corresponds to the contributions to the Josephson current coming from Equation (|47[) . i.e. from inter- 
ference of pair tunneling and the initial population, 

M ^ = WT^) + WT5^) (50) 

and Mgingie to that coming from Equation (|46[) . i.e. from the interference in the single-particle tunneling to the 
excited state 1 14) 

M **&- (u-5 M )(u + s 13 y (51) 

From equation (|49j) we obtain the relation given in the article for the Josephson current 

1 1 
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sm[(5 13 + 6 24 )t + ip]. (52) 
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Coupling dependence of the Josephson frequency 

The exact numerical solution of the four-state model is very much feasible, and supplements well the pcrturba- 
tive analytical treatment above. In Figure [5] we show that for small couplings the perturbative treatment agrees 
with exact numerics as one would expect. We may then study numerically the higher order effects in the coupling 
strength. The dependence of the Josephson frequency on the coupling strength is a particularly interesting question. 
In Figure [7] we plot the value of the Josephson frequency as a function of the RF coupling. As mentioned in the 
article this coupling dependence of the Josephson frequency may impose a limit to the accuracy of the Josephson 
voltage-frequency relation. 

Numerically calculated eigenstates 

In general, solving analytically the eigenstates of the four state system of Figure [5] produces extremely cumber- 
some formulas which do not give much insight to the problem. However, to illustrate our arguments, we give here 
numerically calculated eigenstates in certain selected cases. Since we are now solving the eigenstates of a multiple 
state system where the states are coupled and have energy differences (detunings) between them, this is similar to 
dressed state descriptions of some quantum optics systems pQ. 

Zero detunings 

Let us consider a system with U/J = —5.0. In case of zero detunings S\ 3 /J = 624/ J = but having finite 
couplings, here ^13/ J = ^24/ J = 1, the dressed states are the following 

|7) = -0.6672 |12) + 0.2341 114> - 0.2341 123) - 0.6672 |34) , (53) 

|7/> = 0.7071 ]12> — 0.7071 [34) , (54) 

\III) = 0.7071 [14) + 0.7071 123) , (55) 

\IV) = -0.2341 112) + 0.6672 [14) - 0.6672 |23) - 0.2341 [34) , (56) 
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Figure 6: The Josephson critical currents |J^| (dots) and |i^| (crosses) as a function of S24/J, when 813/ J = 2.0, 



given by the exact numerical calculation for the four-state system (here U/J 
show the results of the perturbative calculation. J is the unit of energy. 



5.0 and O/J = 0.1). The solid lines 
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Figure 7: Position of the Josephson peak (i.e. the Josephson oscillation frequency) as a function of the RF coupling, 
with ft 13 = n 2i = O, 5vi /J = 5 2 i/J = 1.0 and U/J = 10.0. 



with the eigenencrgics 

Ei = -5.7016, E n = -5.0, E m = 0.0, 7? /v = 0.7016. (57) 

In general, if the initial state of the system is one of eigenstatcs of the Hamiltonian (a dressed state), there will be no 
dynamics. Here the state \II) is what one could consider in our small system as an analogue of the BCS state for two 
supcrfluids: a superposition of pairs. Having that state as the initial state, no dynamics takes place; this corresponds 
to the absence of Josephson oscillations for zero voltage (zero detuning) in case of identical superconductors. But if 
the phase difference between the superconductors, or their other properties, would be such that the initial state is 
no longer identical to one of the eigenstatcs (e.g. \II)), then there will be dynamics. That corresponds to the DC 
Josephson effect, caused by a phase difference between two superconductors instead of a voltage. This is just the 
standard two-mode description of Josephson physics found in many text-books (the two modes correspond to the 
pairs (superconductors) on both sides of the junction, here to the states |12) and 1 34) ) . 

Symmetric detunings 

Now let us look at the effect of the detunings when they are symmetric. We take here the values S13/J = 624/ J = 1. 
Then the eigenstates are 

\I) = 0.9581 112) - 0.1737 |14) + 0.1737 |23) + 0.1470 |34) , 
\II) = -0.2072 |12) - 0.17 1 14) + 0.17 |23) + 0.9483 [34) , 
\III) = 0.7071 [14) + 0.7071 123) , 

\IV) = 0.1976 [12) + 0.664 [14) - 0.664 |23) + 0.2813 |34) , 

with the eigenencrgics 

Ei = -6.3626, E u = -4.3586, Em = 0.0, E IV = 0.7212. (62) 

Now one can see that the states 1 12) and |34) do not appear any more symmetrically in the eigenstates. This leads 
to the Josephson oscillations for finite detuning (finite voltage), if we have an initial state that contains only 1 12) 
and 1 34), as would be the analogue for superconductors. 

Note that the eigenstates contain unpaired states. These can then contribute to the dynamics. In this sense, the 
four state system description goes beyond the simple two-mode text-book description of the Josephson effect. Of 
course, these unpaired contributions are implicitly involved in more elaborate descriptions of Josephson physics, as 
is for instance evident from the fact that our self-consistent linear response calculation could reveal the asymmetry 
of the critical currents in the asymmetric detunings case. 

However, in previous literature, the unpaired contributions have not been paid attention to, because, as is clear 
from the symmetric structure of the above eigenstates with respect to states 1 14) and |23), the contributions of 
these states will be equal in magnitude in the dynamics. Therefore it appears as if only pairs are tunneling and no 
single particles exist during the dynamics: but in fact, the single particles exist, it just happens that the expectation 
values of single particles in states 1 14) and 1 23) at any given time of the dynamics are the same. Therefore it seems 
like |1) and |2) arc tunneling together. Another reason why the single particle contributions have not been notified 
earlier is that it is known that the standard single particle currents require detunings (voltages) above twice the gap; 
how could they then participate in oscillations at the below gap Josephson frequency? In this manuscript, we have 
shown that it is the interference term of these single particle currents (i.e. a beating of the standard single particle 
contributions) that contributes to the Josephson oscillations. 



(58) 
(59) 
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Asymmetric detunings 

Finally, let us choose different detunings with 613/ J = 1 and S24/J = 2 i.e. we consider the "different voltages 
spin up and spin down electrons" case. Then the eigenstates are 

\I) = 0.9704 |12) - 0.1684 [14) + 0.1454 |23) + 0.0944 |34) , (63) 

1 77) = -0.1582 |12) - 0.2291 114) + 0.1774 |23) + 0.944 |34) , (64) 

|777~) = 0.0884 |12) + 0.9138 [14) + 0.3586 |23) + 0.1692 |34) , (65) 

\IV) = -0.1599 |12) - 0.2902 |14) + 0.9049 |23) - 0.2672 |34) , (66) 



with the eigenencrgics 



Ei = -6.8234, Eu = -3.9305, Em = -0.2182, E TV = 0.9721. (67) 



The observations of the dynamics in the previous case of symmetric detunings apply here as well, with the important 
difference that now the single particle contributions 1 14) and [23) do not appear symmetrically in the eigenstates. It 
is then understandable that asymmetry with respect to them appears also in the dynamics: this is the asymmetry in 
critical currents that we have predicted and explained. Note that by finding and explaining this asymmetry, we could 
provide insight also to the standard symmetric case as explained in the case of symmetric detunings: the Josephson 
effect consists of interference terms of pair and single particle tunneling processes, in particular, the interference of 
standard single particle currents contributes to the Josephson effect also in the symmetric case, although it appears 
as pair tunneling due to the symmetry. In the asymmetric case, the single particle interference term then clearly 
manifests itself, and the novel intuitive understanding of the Josephson effect proposed by us, namely interferences 
in pair and single particle tunnelings, can be put under a direct experimental test. 
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We propose that with ultracold Fermi gases one can realize a spin-asymmetric Josephson effect 
in which the two spin components of a Cooper pair are driven asymmetrically — corresponding 
to driving a Josephson junction of two superconductors with different voltages Vf and Vj, for spin 
up and down electrons, respectively. We predict that the spin up and down components oscillate 
at the same frequency but with different amplitudes. Furthermore our results reveal that the 
standard interpretation of the Josephson supercurrent in terms of coherent bosonic pair tunneling 
is insufficient. We provide an intuitive interpretation of the Josephson supercurrent as interference 
in Rabi oscillations of pairs and single particles, the latter causing the asymmetry. 
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When a coherent many-body system is partitioned into 
two sub-systems, the dynamics of macroscopic observ- 
ables such as relative number of particles and relative 
phase is called the Josephson effect 0, The exter- 
nal Josephson effect has been realized in superconduct- 
ing junctions superfiuid 3 He [B| and 4 He @, and in 
Bose-Einstein condensates (BEC) of alkali atomic gases 
in double- well traps 0, H|. The internal Josephson ef- 
fect has been demonstrated in 3 He Q and is expected 
to occur in spin BECs d, [l(]l - fl2| . Also in the context of 
ultracold Fermi gases [l3| the possibility of the Joseph- 
son effect has recently received theoretical interest 14- 



that the spin-asymmetric effect can be realized in a SIS- 
junction of two materials with different Zeeman splittings 



171 ]. In this letter we show that partitioning a system 



of Cooper-paired fermions so that the two components 
of the pair experience different potentials (this is what 
we mean by "spin-asymmetric" here) leads to a novel ef- 
fect, namely different-amplitude but phase-synchronized 
number-oscillations of the components. Although the mi- 
croscopic description of the Josephson effect is based on 
single particle tunneling, the standard interpretation of 
the Josephson supercurrent is given in terms of coherent 
tunneling of bosons or Cooper pairs [l8j |. Importantly, 
our results show that such an interpretation is insuffi- 
cient. We provide a clear, intuitive explanation of the 
predicted spin-asymmetric Josephson effect and a new 
understanding of the Josephson supercurrent as a pro- 
cess where not only pairs but also the spin-components 
separately contribute via interference. 

We propose that the spin-asymmetric Josephson effect 
can be realized in a four-component Fermi gas in which 
two superfluids are coupled by radio-frequency (rf) fields 
(Fig. la), as in rf spectroscopy [ill]. The setup is moti- 
vated by the recent realization of three-component Fermi 
gases |20j. Another possible, perhaps experimentally 
simpler, realization is a superfiuid two-component ultra- 
cold Fermi gas (Fig. lb) in a (spin)component-dependent 
double-well potential. The theoretical descriptions of the 
systems of Fig. la and Fig. lb are identical (in this let- 
ter we use the notation of the former). We also suggest 



21] 



The setup of Fig. la corresponds to a 

many-body Hamiltonian H = Hq + H r f, 
where H Q = /dr^V ^t(r)(-|i - /i;)V>i( r ) + 
3 / d r UiM (r)Vj 0)^ (r)^ (r). Here ^(r) 
and ^J(r) are the fermionic field operators for the 
internal state i and fii is the chemical potential (we 
assume /i^ = /i), and Uij give the interaction strengths in 
the s-wave contact potential approximation. We assume 
that U12 and U34 lead to pairing and set h = I. 

In the rotating wave approximation [22] ] the tunnel- 
ing coupling between states \i) and \j) is given by the 
Rabi frequency fijj with the detuning Sij = Vij — 
ujij. Here Uij is the frequency of the field and ujij 
is the resonance frequency of the hyperfine transition. 
The effect of the electromagnetic field on the system 
is then described by H r f = ^ J dr(ipl(r)i(j 1 (r) — 

^(r)^(r)) + 5 -f f dr(^(r)^(r) - VlW^W) + 
fl 13 / drip{(r)ijj 3 (r) + h.c. + fl 24 j dr V4(r)^ 4 (r) + h.c. 
In analogy to the usual Josephson junctions, the states 
|1) and 1 2) correspond to spin up and down electrons in 
the left side superconductor, |3) and |4) on the right. The 
detunings £13 and S24 play the role of the voltage. 

The essential new feature in atomic gases is the possi- 
bility to set <5i3 ^ S24, which in the case of the Josephson 
junction corresponds to spin-dependent voltages. Note 
that this is different from superconductor-ferromagnet- 
superconductor (SFS) structures [23| in which the spin- 
active barrier coupling plays the crucial role. Though in 
our case also the couplings could be different, only the 
spin-asymmetric potential is relevant for our predictions. 
Moreover, one might consider ferromagnetic supercon- 
ductors [24j in a junction as a related system but those 
materials have most likely an exotic ground state which 
does not fit our description. 

We now determine the transition rates (i.e. particle 
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currents) between states |1) and |3), /13(f) = (JVi), and 
between states |2) and |4), /24(f) = {-N2). We calculate a 
self-consistent linear response with respect to H r f (valid 
when the number of transferred particles is small com- 
pared to the total particle number) for the system of Fig. 
1 with the aid of the Kubo formula and the Kadanoff- 
Baym method [25;]. We obtain the currents 



/13(f) = I13 + l! L3 shi[(<5 13 + 5 2i )t + ip], 
/24(f) = I24 + /£ sin[(<5i3 + Mf + <P]- 



(1) 



Here I s is the standard single particle (quasiparticle) cur- 
rent that occurs only for detunings Sij above the excita- 
tion gap 2A. The initial phase of the Josephson current 
is ip. The critical Josephson currents I c become, in a 
spatially homogeneous case and in the BCS description, 

ig =2 \n 13 n 2i u^( P = 0, 5 24 + iv + )\ , (2) 
ig =2 |r! 13 r! 2 4iMp = o,Ji 3 + l? 7+)| , (3) 

where n>(p,w) = ^ J2^ x ^"12 (q, x)FtS?L ~ P) X - w). 
Here V is the volume and (3 — l/(fcsT) (T is temperature 
and ks the Boltzmann constant). T\2 = (.7-34 = .F43) 
is the anomalous mean field Matsubara Green's function 
for the superfiuid of components |1) and 12) (|3) and |4)). 
For details see supplementary material [26j] . 

The striking result is that the critical current /f^ in 
equation $2$ depends only on 824, and similarly ig in 
equation ^ only on #13 (this is not limited to the BCS 
regime; it remains true whenever pairing correlations ex- 
ist). By choosing different detunings (5 13 and 624, one 
can observe a spin-asymmetric Josephson effect in which 
the currents in the two tunneling channels are different 
in amplitude, but oscillate at the same frequency. More- 
over, the results predict a tunable DC Josephson effect: 
by choosing the detunings so that S13 + S24 = 0, the phase 
factor in equation ([1]) is constant but the critical current 
can still be tuned. The conclusions hold for experimen- 
tally realistic parameters and also if cross interactions 
(see Fig. la) and finite temperature are included into our 
analysis, as shown in the supplementary material (2(|. 
For typical parameters and taking e.g. Si 3 /E F = 0.4 
and 824/Ep = 0.5, with Ep the Fermi energy, one ob- 
tains a considerable asymmetry of /13//24 = 1.14, see 
[26l | for details. By performing the self-consistent cal- 
culation we have removed the ambiguity of whether our 
previous suggestion of the critical current asymmetry 14] 
was due to a simple linear response approach following 
the Ambegaokar-Baratoff treatment [27[ . 

Our results CD)-© arc in obvious contradiction with 
the standard interpretation of the Josephson super- 
current in terms of coherent pair tunneling exclud- 
ing any difference in the Josephson currents of differ- 
ent spin-components. In what follows we provide an 
explanation for the spin-asymmetric Josephson effect 
which not only resolves this paradox but also opens 




FIG. 1: Spin-asymmetric Josephson effect setup, (a) The 
four-component Fermi gas. Particles in internal (e.g. hyper- 
fine) states |1) - 1 4) coexist spatially. The components |1) 
and 1 2} as well as |3) and |4) form Cooper pairs due to the in- 
teractions U12 and (734. The cross-interactions U13, U14, U23, 
U24, are assumed weak, thus not leading to pairing. The tun- 
neling is driven by applying rf fields to couple the states. The 
Rabi couplings Qij correspond to the weak tunneling link be- 
tween the two superconductors in a Josephson junction. The 



detunings Si 



uiij play the role of the voltage, eV. 



Remarkably, one can choose 613 7^ 824 and create the ana- 
logue of a spin-dependent voltage in a Josephson junction. 
The states |1) and |3) (|2) and |4}) must be states of the same 
atom to allow the rf coupling, but |1) and |2) can be e.g. 6 Li 
and 40 K. (b) Cooper pairs of a two-component Fermi gas in 
a spin dependent double- well potential. This setup is concep- 
tually equivalent to the system of (a) albeit the absence of 
cross-interactions. 



up a new point of view to the conventional Joseph- 
son effect. Let us begin by writing down the ini- 
tial state of the two superfluids as a product state 
of two BCS states: |*) = |BCS} 12 cg> |BCS) 34 = 
Uk K |0> fc + v k 1 12},) n* («*' \0)v + «*' |34} fe ,) where 
\ij)k = c k,i c -k,j Next, we single out one Cooper 

pair from each superfiuid. Since the rf coupling between 
1-3 and 2-4 conserves momentum we focus on states with 
k = k' (the momentum conservation can be relaxed and 
our conclusions still hold): 

K |0} fe + v k |12} fe )K |0) fe + v k |34) fe ) = u\ |0) fc |0} fc 
+ v\ |12) fc |34) fc + u k v k \12) k \0) k + u k v k |0) fe |34) fc . (4) 

The empty state |0) |0) cannot contribute to the current 
and neither can 1 12} 1 34) since it is Pauli blocked. There- 
fore, the Josephson physics arises from the u k Vk 1 1 2) |0) + 
UkVk |0) 1 34) superposition. 

We then ask whether the essential features of our 
results can be explained by considering the dynamics 
of a single Cooper pair, initially in the above super- 
position state characteristic for the BCS state with a 
macroscopic phase. In additition to the paired states 
12) |0) = |12) and |0) 1 34) = |34) the states |1) |4) = |14) 
and 1 2) |3) = 1 23} are required to have a closed subsys- 
tem with respect to the tunneling coupling, see Fig. 2. 
These broken pair states are analogous to single particle 
excitations within the BCS formalism. 

We now solve the time evolution of this system per- 
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turbatively in the couplings Oy. We consider an initial 
state in the general superposition form ct\ 1 12) + (3u |34) 
as suggested by equation (@} • The total current becomes 



<0(*)| JVi |0(*)> 



(Mt)\fii\Mt)) 



■\Pii\ 2 (Mt)\Ni\Mt))- 
•ofAi (<M*)l^il<M*)>, 



(5) 



where |<fa(*)> = exp(-ifft) 1 12) and |0n(t)) = 
cxp(-iHt) 1 34) are calculated to second order in Qy. For 
details and for a complementary discussion in terms of ex- 
act numerical eigenstates (dressed states) see [26}. Here 
the first two terms contain only single particle Rabi pro- 
cesses, which correspond to the standard single particle 
(quasiparticle) currents in a Josephson junction, in 
equation ([T]). Only the last two terms in equation ([5]), 
with (XjPij = \ai/3u\e l,p , contribute to the Josephson cur- 
rent. Thus, the Josephson effect originates from the in- 
terference part of the Rabi oscillations in the |12), |34), 
1 14), 1 23) state space. 

Isolating the terms oscillating at the Josephson fre- 
quency in {4>(t) \ Ni \(j)(t)), we get the Josephson current 



I13 = sin[((5 13 + 5 2 i)t + ip], 



lg(S 2 4)=2il 13 n 2i \a l p lI \ 



1 



1 



U + 5. 



2 I 
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(6) 



Hereby, we obtain the same qualitative result as our lin- 
ear response calculation gave in equation (fTJ : the ampli- 
tude of Ii 3 depends only on the detuning 824 ■ Now we can 
identify the source of the dependence of lf 3 on 624- The 
current lf 3 derives from the population of species |1), i.e. 
both the population of state |12) and |14). The result of 
equation © is the sum of these two contributions. 

The contribution from state 1 12) is the result of tun- 
neling between the paired states 1 1 2) and |34), via the 
intermediate state 1 14) (or 1 23) ) as shown in Fig. 3a. 
The term is symmetric in $13 and 624 and proportional 
to M vn ir = -rrrs h n h ■ The contribution is of sec- 



ond order because it originates from a zeroth order and a 
second order process: the population due to pair tunnel- 
ing from 1 34) to 1 12) (second order) is indistinguishable 
from the initial (zeroth order) population in state 1 12) . 
Pair tunneling contributions that do not originate from 
interference also exist but they are of fourth order in ft. 

The contribution to /f^ from state 1 14) arises from in- 
terference of broken pairs (single particles) as depicted in 
Fig. 3b. This contribution is responsible for the asym- 
metry as it is proportional to M sing i e = v } S2i - u+ &13 ■ 
Physically, this term arises because we cannot distinguish 
between broken pairs from the states 1 12) and |34). The 
two paths leading to the state 1 14) are not symmetric 
with respect to 613 and (524, see Figs. 2 and 3. This 
gives important new insight also to the Josephson effect 
in general: also single particle interferences contribute to 
the Josephson current if, sin(<5i 3 + S24) of equation (fT|). 
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FIG. 2: Energy level diagram of the four-state model. The 
basis states |12), |34), |14), 1 23) form a closed subspace under 
the Rabi oscillation processes caused by the rf coupling. The 
paired states 1 12) and 1 34) have an energy lower by U than 
the states of broken pairs 1 14) and 1 23) . 
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FIG. 3: Interference processes leading to the Josephson su- 
percurrent. The initial state is \<f>) = (|12) + |34))/V2. (a) The 
second order transition from 1 34) through 1 14) to 1 1 2) creates 
an interference with the initial (zeroth order) population in 
state 1 12). The resulting contribution is of the order 013^24- 
Also the same process via the state 1 23) contributes, (b) The 
first order transitions to the excited state 1 14) from 1 12) and 
1 34) create an interference proportional to f2i3Sl24- 



In the standard symmetric case 643 — S24 = S, we have 
M pair = 2/(U + S) and M smgle - 26/ {U 2 - 5 2 ). At 
the DC Josephson limit, 5 — > 0, the pair transfer M pa i r 
dominates, which is intuitively appealing. Closer to the 
Riedel peak the excited state (single particle) interfer- 
ence M s i n gi e becomes equally important. To illustrate 
further, in a typical Al/AlO^/Al junction at the volt- 
age of 0.015 mV the single particle interference accounts 
for 1.7 % of the AC Josephson current. Note that our 
" single particle interference term" is not the cosine-term 
(also called "quasiparticle interference term" [2^]) of the 
Josephson effect. The cosine-term involves real single 
particle transitions and exists only for voltages above 2A 
at zero temperature. In contrast, our single particle in- 
terference term is inherent in the supercurrent and corre- 
sponds to virtual transitions. Similarly, our findings are 
different from various combined effects of single particle 
currents and supercurrents in small Josephson junctions 
[29I ] . where again the single particle transitions are real, 
not virtual. Note also that we do not consider any inter- 
actions between the Cooper pairs (analogue of charging 
effects) nor the effect of the environment. 

We emphasise the interference nature of the Joseph- 
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son effect: it requires lack of which- way information 
on the tunneling path of the particle. Above, the 
superposition u k v k (\12) |0) + |0) 34}) was due to the 
uncertainty of the particle number in the BCS state. It 
is interestin g to contrast this with the number-projected 
BCS-state |30j or the Fock-state 31], \^) Ncons = 
(u k v k , |0} fe \\2) k , + u k ,v k |12) fe |0) fe ,)( TO |0} fc |34} fe , + 
u k 'V k \0) k , |34} fc ) (for two ^-states.) Here the 
part relevant for Josephson physics would be 
u k v k 'U k 'V k |12) fc , \M) k + u k v k >u k >v k \l2) k |34} fc , . Now, the 
entanglement of |34) fe , with 1 1 2 ) ^ allows to determine, 
whether a particle in state |1) belongs to an initial 1-2 
pair or to a tunnelled 3-4 pair. The indistinguishability 
is lost and there is no Josephson effect. Note that 
the system of Fig. la should be cooled with the rf 
couplings on to realize the necessary uncertainty in 
particle number. One can also ask whether, in the case 
of the separated 1-2 and 3-4 condensates (Fock states), 
the measurement process itself is sufficient to generate 
the relative particle number uncertainty, in analogy to 
the famous problem of interference between two BECs 

& 

Importantly, based on the picture of interfering Rabi 
processes, one would anticipate corrections from the cou- 
pling strength fly to the Josephson frequency 5\z + ^24 
given by linear response. We have observed such a cor- 
rection in the numerical simulations, see supplementary 
material [2q |. We expect this correction to vanish in 
the thermodynamic limit, but it could be experimentally 
tested in ultracold gases and in superconducting grains. 

In summary, we have predicted a spin-asymmetric 
Josephson effect which could be observed in ultracold 
Fermi gases and solid state systems and is likely to be 
relevant for spintronics. We provide a microscopic de- 
scription of the Josephson effect as interference in Rabi 
oscillations of pairs and single particles, where the latter 
cause the predicted asymmetry. In particular, our find- 
ing that the interference part of single particle currents 
contributes to the Josephson supercurrcnt, in addition to 
the interference in pair tunneling, is fundamental. The 
single particle interference contribution is present already 
in the standard (symmetric) case, and becomes manifest 
and experimentally verifiable in the asymmetric one. 
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